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Abstract 

The Donnachie-Landshoff fit of total cross sections has now be- 
come a standard reference for models of total, elastic and diffrac- 
tive cross-sections. Adopting their philosophy that simple-pole 
exchanges should account for all data to present energies, we as- 
sess the uncertainties on their fits. Our best estimate for the 
pomeron intercept is 1.096io]oo9) but several models have a good 
X 2 for intercepts in the range [1.07,1.11]. 
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In Regge theory, the energy dependence of total and elastic cross sections 
is implied by the analytic structure of the hadronic amplitude. The simplest 
singularities are poles in the J-plane, which correspond to the existence of 
bound states of mass M and spin J with J = a(M 2 ). The sum of the 
exchanges of all the particles on a trajectory a(M 2 ) leads to a hadronic 
amplitude which behaves, for s —>■ oo and t = M 2 < finite, as: 



The last term, which arises from s to u crossing, is usually included as a Regge 
signature factor. Its sign depends on the charge parity of the exchange. 

Meson trajectories are observed to be linear, as can be understood in a 
string model: 



with a the intercept, and a' the slope of the trajectory. In particular, total 
cross sections are related to ImA(s, t = M 2 = 0) by the optical theorem 
and involve only the intercept of the trajectory. Many such trajectories are 
present, but at high s only the highest-spin trajectories will contribute: the 
others are suppressed by a power of s. Hence, from the observed meson 
spectrum, one expects total cross section to fall like l/y/s, as the leading 
meson trajectories are those of the p, u, a and / mesons, with a ~ 0.5. 

The fact that total cross sections rise at high energy forces one to in- 
troduce one further ingredient: the pomeron, which in its simplest version 
is another Regge trajectory, with an intercept slightly greater than 1. This 
object has been hypothesized a long time ago Jl|], and presumably arises from 
the gauge sector of QCD. Its intercept and slope are then fundamental num- 
bers characterizing the pure gauge sector of Yang- Mills SU(3). They control 
both the large-s and the small- x limits of hadronic cross sections. 

There has been a renewed interest in this object as the range of s and 1/x 
has been expanding, at HERA and at the Tevatron, where e.g. rapidity gaps 
in deep inelastic scattering have been observed [0. The situation at HERA 
is particularly confusing, as the pomeron intercept seems to vary with the 
scale of the problem, and with the masses of external particles: whereas 
it seems to be of the order of 1.1 in the total 7p cross section, it increases 
slightly in quasi-elastic production of p mesons, reaches values of the order of 
1.25 in J/ip photoproduction, and increases further in high-Q 2 deep-inelastic 
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scattering. At the same time, the gap cross section remain compatible with 
a low value of the intercept, as measured in the pomeron flux. 

The only object which is well defined operationally is the soft intercept, 
i.e. the one which controls total cross sections at high-energy. How well that 
number is known is thus a central issue for HERA, and will no doubt have 
a bearing on the extrapolation of soft cross sections to LHC energies. Hence 
we propose in this letter to re-evaluate the simple-pole parametrisation and 
to estimate the errors on the various parameters. 

The information on the soft intercept comes from total cross sections. As 
the pomeron controls the high-energy behaviour of the cross section, it will be 
most sensitive to the highest energy (i.e. pp and pp) data. There is one caveat 
to this, as one expects the model to break down at low energy (because of 
the presence of sub-leading trajectories) and at high-energy (because multiple 
pomeron exchanges must unitarise the cross sections to make them agree with 
the Froissart-Martin bound). Hence one must not only determine the best 
parameters, but also the range of validity of the model. One might think 
of resorting to a \ 2 test i n order to do this. This is not possible directly 
because of the low quality of the data. Indeed, the largest compilation of 
data, made available by the Particle Data Group || contains points which 
are inconsistent at the 16<t level! [Note that this point is not shown any longer 
in the 1996 version of the particle data book, although it is still present in 
the compilation of data]. 

To illustrate these problems and their proposed solution, we first examine 
the Donnachie-Landshoff (DL) model in its original form ||. This model fits 
the pp and pp cross sections using a minimal number of trajectories. The 
higher meson trajectories, a/ f (C=+l) and p/uj (C = —1), are assumed to 
be degenerate, and their contributions are added to a pomeron term in the 
amplitude. The data are fitted for y/s > 10 GeV, as lower trajectories would 
then contribute less than 1%, which is less than the errors on the data. The 
result of their fit [[J is a pomeron intercept of 1.0808, for which they did not 
quote a x 2 or error bars, but simply mentioned that the y 2 was very flat near 
the minimum. 

We show in the first column of Table 1 our results for such a fit. We use 
the usual definition 
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with di ± ei the measured pp or pp total cross section at energy y/sl, and 

ImA(s, 0) = C_s a ™ + C + s am + C P s ap (4) 



where C_ flips sign when going from pp to pp. We shall use the notation 
«p = 1 + e in the following. 



parameter 


all data 


filtered data (2a) 


filtered data (la) 


x 2 

X 2 per d.o.f. 


410.8 
3.16 


80.3 
0.62 


32.4 
0.25 


pomeron intercept-1 
pomeron coupling (mb) 


0912 +u - uuvv 
u.uyiz_ .0 070 

19.3ti7 


0.0887™ 
19.8l{- 7 


0.0863™ 
20.4+^ 


p/u/a/f intercept-1 
pjijj coupling C-(pp) (mb) 
a/ f coupling (mb) 


n QOO+U.U65 

u.ooz_ 071 


n Q7Q+U.U67 
u -°' J -0.073 

-13-3^1 
62111 


-u.dy»_ 090 

-15 3+ 5 - 7 

fi7 +25 
°' -16 



Table 1: Simple pole fits to total pp and pp cross sections, assuming 
degenerate C = +1 and C = — 1 exchanges. 



The value of the \ 2 = 410 for 135 data points and 5 parameters is totally 
unacceptable, as it corresponds to a confidence level (CL) of 2 x 10~ 36 ! This 
CL gives the probability that the model with fixed parameters could generate 
the data through a random fluctuation. Note that all other fits have so far 
suffered from the same problem || . There are two possible outcomes to such 
a high value of the \ 2 '- either the model is to be rejected, or some of the 
data are wrong. As we already mentioned, there are a few obviously wrong 
points within the data. Hence before reaching conclusions about the model, 
one needs to tackle the issue of eliminating those points. 

We propose here a reasonable criterion, which will give us the order of 
magnitude of the uncertainties. The best selection criterion would based on 
physics arguments to decide which experiments got the wrong results. Unfor- 
tunately, besides experimental questions, this would involve some personal 
bias, and prove infeasible for old (ISR) data, where most of the incompati- 
bilities lie. Here, we propose the following approach, which is independent of 
any underlying theoretical model. Most data are bunched in several energy 
intervals. We ask that a given data point should not deviate by more than 
1 or 2cr from the average of all data in a bin of ±1 GeV centered around 
it (note that filtering at the la level rejects both the E710 and the CDF 
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measurements, whereas both are kept when filtering at the 2a level ). All 
the data which do not fulfill this criterion are to be rejected. This eliminates 
some of the pathological points in a fit-independent manner, and brings in a 
dramatic improvement in the x 2 /dof., which falls down to less than 1 for the 
best fits. Note that this procedure is not unlike the one followed by UA4/2 in 
U . The number of standard deviations, as well as the width of the intervals 
considered are of course somewhat arbitrary. Because the data are concen- 
trated at certain energies, the width does not matter too much, but one could 
certainly change the filtering to an arbitrary number of standard deviations. 
When considering large numbers of points, it is in fact quite likely that about 
one third will deviate by more than one standard deviation. Hence the la 
filtering is presumably too stringent, although we shall see that the results 
are stable when going from one selection criterion to the other. 

We give in Table 2 the number of points kept when filtering the data, 
the full data sets being available at |http: / /nuclth02.phys.ulg.ac.be/DATA.html . 
Note that in the following we shall also consider an alternative data set 
Its main difference with that of the Particle Data Group is that the statistical 
and systematic errors have been added in quadrature, and that it includes 
the measurements of the p parameter. 



data set 


Vtot ( m b) 


°tot ( m b) 


pPV 


pPP 


P.D.G. 


94 


28 






P.D.G. @ - 2a 


84 


28 






P.D.G. S - la 


65 


20 






Ref. 1 


66 


29 


41 


13 


Ref. - 2a 


60 


28 


38 


13 


Ref. 0- la 


53 


19 


31 


13 



Table 2: The number of points kept after data selection, for 
Vs > 10 GeV. 



If we now go back to the first line of Table 1, we see that the elimination of 
a handful of points drastically changes the x 2 ■, and makes the model entirely 
acceptable. We can now proceed to evaluate the errors on the N parameters 
[§, 10, ITU . Here the data are given, and we are asking what is the effect of 
a random change of parameters. The value of x 2 ~ Xmin * s distributed as a 
X 2 distribution with N parameters. We choose here to quote the intervals 
corresponding to the projection of the constant x 2 hypersurface containing 
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70% of the probability. This corresponds to x 2 



X,r 



+ 6.06 in the DL case 



with 5 parameters. Note that we do not trust the validity of the method 
used by the Particle Data Group in their fit to the DL model, as they simply 
renormalise the x 2 to x 2 /dof. = 1 and let the new \ 2 vary by one unit. 

An important test of the data filtering method is that the central values 
and their errors should not depend too much on the filtering itself. Table 2 
shows that this is indeed the case for our fit. We see that the pomeron 
intercept is determined to be about 1.090, and that it could be as high as 
1.096. Filtering the data does change the value of xLini but does not affect 
its variation around the minimum too much. Hence the stability of the 
parameter values seems a more reasonable criterion than the value of Xmin- 

We can now tackle the question of the energy range of validity of the 
model. The two basic requirements are that the x 2 /dof. be of the order 
of 1, and that the determination of the intercept be stable. We show in 
Fig. 1 the result of varying the energy range. Clearly, the lower trajectories 
seem to matter for \/s rnin < 10 GeV, whereas the upper energy does not 
seem to modify the results (in other words, there is no sign of the onset of 
unitarisation). Hence we adopt \fs min = 10 GeV as the lowest energy at 
which the model is correct. This happens to be the point at which xLin is 
lowest, too. This dependence on the lower energy cut explains why both the 
Particle Data Group M and Bueno and Velasco ]T^] obtain an wrong value 
for the intercept, much lower than ours. 



<=. Xljd.o.f. 




xliJd-o.f.(2a) X 0.1 

" e(2a) £-1- 



10 




(GeV) 



Figure 1: DL intercept-1 as a function of the lower (a) and upper (b) 
energy cuts on the data. The curve shows the x 2 /dof. for data filtered 
at the 2a level. 
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One must wonder if it is possible to get a better determination of the 
soft pomeron intercept by using more data. After performing their fit on pp 
and pp cross section, DL extended it to include all measured hadronic total 
cross sections. They found that this did not affect the value of the intercept 
very much. This is due to the fact that the other total cross sections have 
large errors, and that new parameters are introduced (namely the couplings 
of the Regge trajectories) each time one considers a new type of cross section. 
The Particle Data Group (PDG) || obtained very narrow determinations of 
the pomeron intercept from the other hadronic reactions. We believe that 
their conclusions are wrong, and illustrate this in the case of the 7r ± p total 
cross sections, for which they use \fs min ~ 4 GeV and obtain an intercept of 
1.079 ±0.003. 

M (b) 




V^™ (GeV) (GeV) 

Figure 2: Pomeron intercept from 7rp data: (a) shows the values of the 
intercept as the lower energy cut on the data is changed; (b) shows 
our best fit together with that of the Particle Data Group, together 
with the set of data filtered at the 2a level. 
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We show in Fig. 2 our results for such a fit: for data filtered at the 2a 
level (92 points), and for y/s min = 4 GeV, we obtain «o = 1.115lo;o23- There 
is no significant change if we modify the lower energy cutoff on the data. We 
show in Fig. (2.b) our best fit together with that of the Particle Data Group. 
Although according to their estimate our central value for the intercept is 10 
standard deviations from theirs, we see that the two fits are indistinguishable. 
Hence we believe that both their standard deviations and their central values 
are wrong. These conclusions are not affected by the use of the full data set 
instead of the one filtered at the 2a level. The above value of the intercept in 
fact gives a slightly smaller x 2 /dof than the one quoted in the Particle Data 
Book. Note that this intercept is consistent with the one we got from the 
analysis of pp and pp total cross sections. The conclusion from this exercise is 
that the errors from the low-energy hadronic data are large, especially if we 
use a low-energy cut-off of the order of 10. Hence we want to limit ourselves 
to pp and pp amplitudes at t = 0. 



parameter 


all data 


filtered data (2a) 


filtered data (la) 


x 2 


561.3 


168.3 


94.9 


X 2 per d.o.f. 


3.28 


1.07 




pomeron inter cept-1 
pomeron coupling (mb) 


0.0840 ± 0.0050 
20.8 ± 1.1 


U.U8±7_ 00053 
21.4 ±1.1 


0.0804™ 
21.8 ±1.3 


p/u/a/f intercept-1 
p/u coupling C-(pp) (mb) 
a/f coupling (mb) 


-0.408™ 
-14.0±il 
67.01™ 


-0 421+ au34 
u.4zi_ 036 

-i6.5±t:l 

66.6^;| 


-0 431+ U - U3V 

-i7.7±i: 9 7 

67.61?;° 



Table 3: Simple pole fit to total pp and pp cross sections, and to 
the p parameter, assuming degenerate C = +1 and C = — 1 meson 
exchanges. 



One more piece of data can be used however: the knowledge of the in- 
tercept is sufficient to determine the value of the real part of the amplitude, 
using crossing symmetry, and hence the measurements of the p parameter 
provide an extra constraint. We use the data collected in Ref. 0], and obtain 
a somewhat worse fit, shown in Table 3, even when filtering data at the 1 
or 2a level. For the 2a filtering of the data, the confidence level goes from 
99.4% to 36%. We show the curves corresponding to the second column of 
Table 3 in Fig. 3 with dotted lines. Whether one should worry about this, 
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and about the change of central value for the parameters, is a matter of taste. 
Note however that it is this small change of central values, combined with the 
effect of too low an energy cut, that lead Bueno and Velasco |T2|] to conclude 
that simple-pole parametrisations were disfavored. 





degenerate 
non degenerate 



\/? (GeV) 



(GeV) 



Figure 3: Best fits to 2a filtered data. The dotted lines correspond 
to the original DL model given in Eq. (||), whereas the plain ones 
correspond to a model where the degeneracy of the lower trajectories 
is lifted, as in Eq. (|6]). The data points are the PDG data filtered at 
the 2a level. 

Before concluding on the best value of the intercept, we need to examine 
the influence of low energy physics on the determination of the intercept. 
Although the energy cut \fs min eliminates subdeading meson trajectories, 
there is still an ambiguity in the treatment of the leading meson trajectories. 
In fact, we shall now see that a slightly different treatment to that of DL 
leads to a better \ 2 an d to more stable parameters. Indeed, there is no reason 
to assume that the p, lu, f and a trajectories are degenerate. We show the 
data for the meson mass spectrum || in Fig. 4, as well as the separate best 
fits. 
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M 2 (GeV 2 ) 



Figure 4: Best fit to the lower trajectories. 
The trajectories plotted correspond to the following equations: 

otu{t) = (0.3467 ± 0.0043) + (0.9213 ±0.0071)* 

a p (t) = (0.5154 ± 0.0014) + (0.8198 ±0.0014)* 

a a (t) = (0.541 ±0.061) + (0.839 ±0.035)* 

a f (t) = (0.697 ± 0.021) + (0.801 ± 0.012)* (5) 

At this point, it is tempting to use the central values of Eq. ([5]) for the 
lower trajectories. Such a procedure produces a a pomeron intercept of 1.095 
with acceptable \ 2 - However, such a procedure is not satisfactory because 
it excludes the possibility of a mixing of the / and the pomeron, or of the 
simultaneous exchange of several trajectories. 

On the other hand, the data are not constraining enough to determine the 
effective intercepts of the four meson trajectories together with the pomeron 
intercept. We adopt an intermediate approach, which is to assume the ex- 
change of separate + and — trajectories with independent intercepts: 

ImA(s, 0) = C^s a - + C + s a+ + C P s ap (6) 

The resulting numbers are shown in Table 4, and are plotted in Fig. 3 with 
plain lines. The \ 2 is smaller than previously, and the parameters are more 
stable. The bounds on the soft pomeron intercept hardly depend on the 
criterion used to filter the data, and intercepts as large as 1.108 are allowed. 
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parameter 


all data 


mterea data (^crj 


filtered data (1<t) 


X 

X 2 per d.o.f. 


DUD. 4 

2.99 


i in 

0.77 


K*7 C\ 
D 1 .0 

0.47 


pomeron inter cept-1 
pomeron coupling (mb) 


0.0990iffl 

Lt.O_ 20 


09fi4+ u - Uiib 
u.uyo4_ 0091 

18.0+1° 


0.095^ 

i8.2il:i 


p/uj intercept-1 
pjijj coupling C-(pp) (mb) 


_n 4Q4+ U - U5ti 
u -* y4 -0.066 

—94 n +6 - 8 

Z4 - u -10.9 


-0.4981-- 
-26.5+I 2 7 5 


U.OJ.U_ 077 

-28.2;? 5 9 4 


a/ f intercept-1 
a/ f coupling (mb) 


n o 19 +U.U5i 
u.oiz_ 052 

56.8l| 7 


-0.315 ±0.058 
54.9+ 7 -° 


-0.324 ± 0.066 
56.2 +9 7 | 


a toi (1.8 TeV) (mb) 


77.6+^ 


76.S+ 2 ; 9 


76.4+jj-? 


a tot (10 TeV) (mb) 


108.4i^; 7 


106.4+^ 


105.4+^ 


a tot (U TeV) (mb) 


115.8l£ 7 


113.5l 7 | 


112 3 +iU ' 8 



Table 4: Simple pole fit to total pp and cross sections, and to 
the p parameter, with non-degenerate C = +1 and C = — 1 meson 
exchanges. 



parameter 


all data 


filtered data (2a) 


filtered data (la) 


x 2 

X 2 per d.o.f. 


197.9 
1.39 


107.7 
0.82 


56.3 
0.52 


pomeron intercept-1 
pomeron coupling (mb) 


u.uyoo_ 0083 
18.4+ 1 ,; 8 , 


0.0940™ 
18.8+ 2 ; 7 


0.095™ 

18.5_ 2 ; 6 


p/uj intercept-1 
p/uj coupling C-{pp) (mb) 


U.O<JO_ 059 

-SI.et™ 


-U.D18_ 0058 
-28.9+? 8 4 


U.D4U_ Q67 

qo c;+8.8 


a/ f intercept-1 
a// coupling (mb) 


n QQQ+U.U54 
U.OOO_ .055 

58.8_ 6 ; 8 


u.ooo_ 057 
61.5 +9 7 ; 8 7 


U.O4O_ gg 

60.4i 7 ° 9 5 


a tot (1.8 TeV) (mb) 


77.3+^ 


77.2 ±2.6 


77.2+^ 


a tot (10 TeV) (mb) 


106.8ltu 


106.31^ 


106.5 + 7 ^ 


a tot (14 TeV) (mb) 


113.9+^ 


113.21^ 


113.5^ 



Table 5: Simple pole fit to total pp and pp cross sections, and to 
the p parameter, with non-degenerate C = ±1 and C = — 1 meson 
exchanges, and using the alternative data set of Ref. (/]] . 



In order to understand better the treatment of the errors, we give in 
Table 5 the result of a fit to the data of Ref. [0 , where the errors have been 
added in quadrature. We see that the results are very stable, especially those 
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for the pomeron intercept. This is similar to the statement that the filtering 
of the data does not affect the parameters very much. 

At this point, the only additional piece of data might be the direct ob- 
servation of the pomeron, i.e. of a 2 ++ glueball. The WA91 collabora- 



tion has indeed confirmed [13| the WA71 observation of the X(1900) and 
showed that it was a I G J PC = + 2 ++ state, /2(1900). Its mass has been 
observed to be 1918 ± 12 MeV. If we assume that this is the first state on 
the pomeron trajectory, and use the determination a' = 0.250 GeV -2 [j|], we 
obtain ap = 1.0803 ± 0.012. This is the value of the intercept for 1-pomeron 
exchange. The intercepts that we obtained in Tables 1, 3, 4 and 5 from scat- 
tering data cannot be directly compared with this value, as they include the 
effect of multiple exchanges, of pomerons and reggeons. But the values we 
have derived are certainly compatible with the WA91 measurement. Note 
however that the conversion of the glueball mass into a pomeron intercept 
relies heavily on the value of the pomeron slope. An intercept of 1.094 would 
be in perfect agreement with the WA91 observation for a slope a' = 0.246 
GeV -2 . Hence it would be dangerous to mix this piece of information with 
the t-channel information. Our best estimate for the pomeron intercept is 
then: 

a P = 1.0964±™o9i (7) 

based on the 2<j-filtered PDG data, in the non-degenerate case. 

Finally, we can place constraints on physics beyond one-pomeron ex- 
change. The first obvious correction at high energy has to be related to 
unitarisation. Clearly, multiple exchanges must tame the rise of total cross 
sections so that they eventually obey the Froissart-Martin bound. Hence one 
would expect the value of the intercept to be an effective one, which decreases 
as y/s increases. We have already seen in Fig. (l.b) that the data does not 
show any sign of unitarisation up to Tevatron energies. One can confirm this 
by introducing a 2-pomeron exchange term in the amplitude. Although little 
is known about such a contribution, there is general agreement that it must 
be negative, and have an intercept lap — 1. We can then fit the data to a 
form: 

ImA(s, 0) = + C+s a+ + C P (s ap + R s a2P ) (8) 

with a 2 p = 2ap — 1. Using 2cr-filtered data, we then obtain an upper bound 
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on the ratio R of its coupling to that of the pomeron: 



1 2 — pomeron coupling | 
1 — pomeron coupling 



< 4.7% (70%C.L.) 



(9) 



Including such a contribution would bring the best value for the intercept of 
the 1-pomeron exchange term to 1.126^o!os2- 

By following the method of Eq. (^), one can obtain bounds on more exotic 
objects. At the 70% C.L., the ratio of the coupling of an odderon to that of 
a pomeron is smaller than 0.1% (the best odderon intercept would then be 
1.105 and the pomeron intercept become 1.099). This would correspond to 
0.08 mb at the Tevatron. As for the "hard pomeron", there is no trace of 
it in the data. Constraining its intercept to being larger than 1.3 leads to 
an upper bound on the ratio of its coupling to that of the pomeron of 0.9% 
(the soft pomeron intercept then becomes 1.065). This would correspond to 
a maximum hard contribution of 19 mb at the Tevatron. 

In conclusion, we have shown that simple pole fits to total cross sections 
are very successful. We show in Fig. 5 the results obtained in this paper 
together with other estimates present in the literature. 

0.15 
0.14 
0.13 
0.12 
0.11 
0.1 
0.09 
0.08 
0.07 
0.06 



WA91 [13] 



this 



P.D.G. [9] 

BV [12] 



DL [5] 



0.5 1 

X 2 /dof 



1 .5 



Figure 5: Our results for the pomeron intercept, compared with others 
in the literature. The values of the x 2 /dof are indicated for the points 
of this work only. 
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All the points from this work have an acceptable x 2 > an d the main dif- 
ference between them is either the filtering of data or the physics of lower 
trajectories. Hence all these estimates are acceptable, and we feel that inter- 
cepts as high as 1.11, and as low as 1.07, are possible. When comparing with 
other work in the literature, we have explained that the use of a small energy 
cutoff leads to smaller intercepts, and reflects the fact that sub-leading meson 
trajectories are to be included [note that the original DL fit || used the same 
cutoff as ours, but used a different definition of x 2 [Pl l- Our errors are much 



larger than those of other estimates because we fully take into account the 
correlation of the various parameters, and because our statistical analysis of 
the data is much more careful than previous ones. In other words, we quote 
the projection of the hypersurface containing 70% of the probability, rather 
than letting the \ 2 simply vary by one unit. For this problem, we believe 
that this leads to much more reasonable error estimates, as explained above. 

Note that these results depend only on pp and pp data. We have argued 
that little could be learned from other hadronic reactions, given that they 
are measured at low energy. In particular, we want to point out that our 
fit to total cross sections of Fig. 3 is indistinguishable from the DL fit for 
y/s < 300 GeV, hence the parametrisation we propose will fit the total ^p 
cross section, as well as the irp and Kp data. 
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